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Abstract 

To Malthus, rapid human population growth — so evident in 18th Century Europe — was obviously unsustainable. In his 
Essay on the Principle of Population, Malthus cogently argued that environmental and socioeconomic constraints on 
population rise were inevitable. Yet, he penned his essay on the eve of the global census size reaching one billion, as nearly 
two centuries of super-exponential increase were taking off. Introducing a novel extension of J. E. Cohen's hallmark coupled 
difference equation model of human population dynamics and carrying capacity, this article examines just how elastic 
population growth limits may be in response to demographic change. The revised model involves a simple formalization of 
how consumption costs influence carrying capacity elasticity over time. Recognizing that complex social resource-extraction 
networks support ongoing consumption-based investment in family formation and intergenerational resource transfers, it is 
important to consider how consumption has impacted the human environment and demography — especially as global 
population has become very large. Sensitivity analysis of the consumption-cost model's fit to historical population 
estimates, modern census data, and 21 st Century demographic projections supports a critical conclusion. The recent 
population explosion was systemically determined by long-term, distinctly pre-industrial cultural evolution. It is suggested 
that modern globalizing transitions in technology, susceptibility to infectious disease, information flows and accumulation, 
and economic complexity were endogenous products of much earlier biocultural evolution of family formation's 
embeddedness in larger, hierarchically self-organizing cultural systems, which could potentially support high population 
elasticity of carrying capacity. Modern super-exponential population growth cannot be considered separately from long- 
term change in the multi-scalar political economy that connects family formation and intergenerational resource transfers 
to wider institutions and social networks. 
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Introduction 

Malthus published the first edition of his essay on limits to 
human population in 1798 [1]. Since then, diverse stresses — 
caused by political violence and marginalization, poverty, poor 
nutrition, and infectious diseases — seem not to have checked 
global growth in the least [2,3]. Only a widening decline in fertility 
has recently begun to slow global demographic growth (Fig. 1) [4]. 
Humanity's modern population rise has profoundly impacted and 
transformed ecosystems around the world [5-9]. The modern 
human population explosion co-occurred historically with what 
has variously been described as the modern technological, 
economic, human capital, and ideological eras [10-17]. It also 
co-occurred with the demographic and epidemiological transitions 
to lower mortality, initially involving geographically patchy 
variation in fertility, followed by recurrent, broadening birthrate 
declines [2,18,19]. The dynamics of the ongoing, globalizing 
demographic transition are usually modeled and discussed — at 
least implicidy — as if they were distinct, driven by qualitatively 
different factors than those shaping pre-industrial population 
dynamics [12,20,21]. Put in Malthusian terms, did the cultural, 



political, and economic dimensions of modernity allow human 
populations suddenly and temporarily to escape earlier environ- 
mental constraints, which had apparendy remained in place well 
into the 18 th Century? Of course, Malthus would skeptically 
expect the answer to this latter question to be "no". Nonetheless, 
in seeing only indications that 19 th Century population would face 
insurmountable limitations to growth, Malthus would surely be 
surprised that global population has roughly doubled three times 
since 1798 (see Fig. 1). 

Cohen's landmark — and remarkably simple — model of coupled 
dynamic change in population and ecological carrying capacity 
supports the plausible claim that the cultural and environmental 
transformations underlying industrialization and modern trans- 
portation and communication technologies had systemic roots 
stretching many centuries prior to Malthus's initial publication of 
his Essay [22,23]. Cohen's model provides an elegant exploratory 
analytical tool for investigating recent global human population 
change. In this article it is argued that an extension of Cohen's 
original coupled-difference-equation model facilitates analytically 
evaluating two alternative hypotheses: (1) that the demographic 
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Figure 1. Semilog plot of historical estimates from 1-1950 CE (various shapes), recent UN census data for 1955-2012 (red asterisks 
with gray background), and 21 st Century projections based on variation in fertility and mortality rate trends (solid red line 
bounded by upper and lower range blue dotted lines), showing a logistic growth pattern with a remarkably rapid acceleration 
during the 19* and 20 th centuries. Although uncertainty in the 21 st Century UN population projects encompasses continued growth to ca. 17 
billion, as well as imminent decline toward ca. 6-7 billion, it is clear that demographic growth began decelerating over the past 20 years, and that 
deceleration is continuing. Data from refs. [4,23,24]. 
doi:1 0.1 371 /journal.pone.01 05291 .g001 



transition was caused by an abrupt systemic transition — perhaps 
driven by a cascading socio-cultural rupture — in the late 18 th or 
early 19 th Century, or (2) that more ancient, preindustrial initial 
conditions determined long-term patterns in the elasticity of 
carrying capacity. In order to consider these possibilities, this study 
carries out a sensitivity analysis of the original Cohen (OC) model 
and the extended "consumption-cost" (CC) model. The aim is to 
compare modeled population trajectories — determimstically 
shaped by interaction with a social-scale-dependent carrying 
capacity, which is elastic in relation to change in population — 
with demographic estimates that are independently based on 
historical data concerning trends in fertility, mortality, land area 
occupied, and population densities (Table 1; see Fig. 1) [23,24]. 
The new — and also very simple — modification of Cohen's original 
model can help to clarify the limits to positive elasticity in carrying 
capacity relative to human population. 

Background 

Cohen's seminal formalization of human demographic growth 
defines scale-dependent social organization of resource extraction 
and processing as a factor that dynamically couples population 
change with changes in environmental constraints on population 



increase. In the OC model Cohen modifies the classic Verhulst- 
Pearl logistic growth equation to define environmental carrying 
capacity as a historically dependent variable — that is, a function 

K(t), so that =rP(t)[K(f)-P(t)] [22,23]. He accounts for 

change in K(t), beginning with the basic assumption that carrying 
capacity is elastic with respect to population size, so that its history- 
dependent dynamic would follow a rate of change proportional to 

dK(t) dP(t) 

that of change in the population, — = c — . In grappling 

with the carrying capacity concept, Cohen is taking on a 
theoretically implied phenomenon so general that — all too 
often — its formalized abstraction obscures, rather than heuristi- 
cally clarifies, the systemic ecological factors shaping population 
growth patterns [25] . What is new is that Cohen concisely models 
endogenously driven logistic population change, whose trajectory 
may be usefully compared with standard Verhulst-Pearl trajecto- 
ries [26]. Basically, Cohen establishes a standard of comparison for 
evaluating whether human population growth is positively shaped 
by social-network, rather than constrained by exogenous, fixed, 
niche-defined limits. It is especially theoretically relevant for 
studying human biocultural evolution that — in defining 
carrying capacity as elastically responding to population change 
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[22,23] — Cohen elegantly supports a dynamic niche construction 
approach [27] to studying our socially intensive, transfer- 
dependent life-history adaptations. 

Indeed, he further suggests that the coefficient c (see above) may 
be usefully defined in relation to P{t). Culturally structured values, 
preferences, beliefs, and interlinked — often competing — institu- 
tions are initial conditions that set a (theoretically) constant 
threshold size, L = c(t)P(i), at which an additional brain or pair of 
hands no longer offers any economy of scale. Thus, with 

c(t) = , carrying capacity change decelerates when P{t) 

ryl) 

surpasses L. Using discrete difference equation forms, the OC 
model illustrates how a coupled system of change in population 
and carrying capacity yields a surprisingly close centennial and 
millennial-scale fit with historical population estimates and 
modern population data (Fig. 2). 

Indeed, with starting conditions set at 1 CE, the OC model 
seemed — at first glance — to account for the super-exponential 
growth that has occurred over roughly the past 200 years. 
Although the world has absorbed a net gain of roughly 1.5 billion 
people since the Cohen equations were published, the OC model 
further presciently supported current United Nations projections 
that global growth rates will actually level off and approach zero — 
or possibly even decline — in this century [4] . 

Modeling Resource-Transfer Impacts on Carrying 
Capacity 

Despite the explanatory promise of the OC model, it has not 
been subjected to thorough theoretical and analytical scrutiny. 
One key limitation in the OC model actually involves an admitted, 
explicit assumption [22,23]: resource-extraction efficiency and 
costs associated with investment in fertility and intergenerational 
resource transfers simply attenuate as P(t) approaches and then 
surpasses L. In the OC system dynamics, as P(f) increases, change 
in carrying capacity approaches zero, reaching a demographically 
stationary, Malthusian steady state. Consequently, once it is 
gained, carrying capacity cannot be lost in the OC model. 

In order to address this issue, the CC model involves a simple 
redefinition of the variable c — that is, the coefficient of relation- 
ship between change in carrying capacity and change in 
population, or the population elasticity of carrying capacity: 

c «=^- flP W" (!) 



The coefficient a is a constant rate of discount on marginal non- 
linear ecological impacts of consumption — including biotic and 
abiotic resource depletion — involved in family formation, invest- 
ment in somatic maintenance, fertility, and transfers to offspring 
and descendants. Relevant values of a are constrained so that 
0<fl«l, reflecting the expectation that consumption costs will 
impact carrying capacity growth only at larger population sizes. 
The non-linear impacts themselves are modeled by power 
coefficient b. As aggregate resource consumption rates increase, 
resources may get depleted in disproportionately positive relation- 
ship to population. Values of b > 0 should yield realistic results. 
With the redefinition of c(t), the differential equation for change in 
dK(t) 

carrying capacity becomes: — -y- — =r[L — aP(t) \. Consump- 
tion costs are expected to rise non-linearly because more than 
individual food and water needs are required to sustain high 
numbers. Also needed — or demanded — are material resources for 
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extraction, processing, transportation, distribution, and even 
material support for family formation, intergenerational transfer 
and consumption practices. Moreover, at large population sizes 
factional and institutional conflicts emerge over land-use for food 
production, transportation and storage; potable water extraction, 
transportation and storage; and non-food and water resource 
extraction, transportation, consumption, and discard. Thus, at a 
certain high-population threshold, each additional person has a 
larger and larger negative impact on resource acquisition and 
distribution, raising the costs for survival, fertility, and intergen- 
erational transfers, while leading in some cases to overexploitation 
of renewable resources and depletion of finite ones. 

The Cost of Reproductive Success in Society 

A second, important area left unexplored is how the OC model 
is defined relative to the human natural intrinsic rate of increase, r. 
Cohen presents his original model for heuristic — rather than 
analytical — purposes [22]. However, the value of r that he uses for 
illustrative reasons scales to an infinitesimal per capita annual 
reproductive rate of roughly 1.5 x 10~ 9 offspring. The life-history 
strategy of Homo sapiens — with distinctively long maternal 
gestation and juvenile growth periods; typically long great-ape 
lactation periods; and a variable age at last reproduction between 
35-50 years for both sexes [28-30] — does limit r max in a 
population with stable age structure with a roughly equal sex 
ratio. Yet, this natural intrinsic rate of increase may be estimated 



to be roughly 10 2 <r max <\0 1 (measured as a continuous rate of 
offspring production per capita per annum) [26,31-33]. This is 
nearly seven orders of magnitude greater than the r-value yielding 
good OC-model fit to independent population estimates. A 
surprising implication of Cohen's emphasis on endogenous 
cultural system growth in carrying capacity is that "r" takes on 
a new, virtually flipped definition. Cohen's revision of the classic 
Verhulst-Pearl logistic growth function swaps "driver" and 
"destination". The basic Verhulstian (and Malthusian and 
Darwinian) expectation is that a high natural intrinsic rate of 
increase is limited by logistic growth to an ecological limit 
determined exogenously to the population itself [34]. Thus, 
natural increase drives growth, and essentially constant exogenous 
ecological limits determine the steady state level. However, in 
Cohen's formulation, substantial population growth only occurs 
when endogenous cultural niche construction processes — possibly 
in combination with exogenous changes, such as climatic 
amelioration — raise carrying capacity, K(t), sufficiently above 
prevailing population, P(f). Moreover, this can only occur through 
prior cultural evolution of the potential for economies of scale, L. 

In fact, the independent historical estimates suggest that, 
following a long period of demographic stasis, global population 
nearly doubled between ca. 900-1300 CE — an interval that 
closely coincides with the brief climatic interstadial known as the 
Medieval Warm Period [35-37] (Fig. 3). Cohen's model supports 
a more plausible explanation of how marginal increases in 
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Figure 2. Logarithmic plot comparison of population-trajectory fit between the standard Verhulst logistic growth model and 
Cohen's discrete-step coupled difference equation model of human population and carrying capacity growth. The original Cohen 
model shows the population trajectory based on the parameter values and initial population and carrying capacity conditions as in ref. [22], The 
Verhulst trajectory is based on the same population initial condition, ,Pic£ = 2.52 x 10 8 , as for the Cohen model; static carrying capacity, K, is set to 
the Cohen model asymptotic equilibrium value P t =K t = 1.555 x 10 10 . The natural intrinsic rate of growth, r, is conservatively set to 0.01, well below 
estimates of human r max [26]. Historical population estimate and UN census data shown are as in Fig. 1. 
doi:1 0.1 371 /journal.pone.01 05291 .g002 
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ecological productivity and reductions in temperature and 
precipitation extremes supported a disproportionate, supra-mar- 
ginal expansion in human numbers, at such a broad geographic 
scale. The prior, independent political-economic emergence of 
potential for economies of scale, L, in different parts of the 
world — across Eurasia, in Africa, the Pacific Islands, and the 
Americas [38] — would have allowed uninterrupted cumulative 
growth in carrying capacity, K(t), which exhibited a positively 
elastic response to Medieval Warm Period climatic amelioration. 

In Cohen's framework, the intrinsic rate of increase, r, is 
implicitly the small amount of population growth that can be 
achieved when prevailing cultural systems limit the social 
mobilization of family formation, resource extraction, and 
consumption to near-replacement levels, so that K(t)is only 
slightly larger than P(t). We may distinguish this cultural, social- 
network-limited rate of increase, r soc , from the natural maximum 
intrinsic rate of increase, r max , so that: 



In general, the human case involves r soc «r mllx x0.l. Now, 
Cohen labeled the parameter c(t) [see equation (1) above] the 
"Condorcet parameter". The formalization of change in carrying 
capacity as exhibiting positive elasticity with respect to change in 
population is a nod to Condorcet's optimistic notion that human 
ingenuity would inevitably solve problems posed by growing 
populations [22]. Although Cohen did not explicitly recognize or 
derive the substantial change in r's definition that he had wrought, 
the relationship term, m, between r max and the much smaller r soc 
may be dubbed the "Marx parameter". As Karl Marx stated in 
Grundrisse, "The human being is in the most literal sense a zoon 
politikon — not merely a gregarious animal, but an animal which 
can individuate itself only in the midst of society. . . Whenever we 
speak of production, then, what is meant is always production at a 
definite stage of social development — production by social 
individuals" [39]. Marx emphasized the potential of social 
relations of production to define and isolate the individual as an 
Aristotelian "political animal", but the political-economically 
mediated, dramatic dilution of r max — a phenomenon preliminarily 
suggested by Cohen's initial work with the OC model — brings into 
focus the profoundly non-individuating systemic tie between self 
and society. Biocultural evolution may be substantially defined by 
the especially complex interconnections in which the individual is 
linked by the family and the wider social networks through which 
family is defined culturally — in the process of family members 
working to obtain, hold, and consume resources, in order to invest 
in fertility and transfer resources to multi-generational sets of 
descendants. 

Long-term systemic change in the functional dilution of 
then, is arguably an important, yet largely unrecognized property 
of socio-politically complex human cultural systems. As Graeber 
has recently noted — in discussing the practical, cosmological, and 
ideological dimensions of kingship, state violence, ritual violence, 
and social order in pre-industrial, non-literate contexts — the 
definition and management of reproductively potent human 
populations as resources to control, sustain, or exploit has likely 
always been a well-focused political concern for factional and 
individual interests in culturally structured social networks [40] . In 
any demographically sustainable human population, individuals 
will face a continuous cost for acquiring and holding resources for 
family formation, investing in fertility, and providing transfers to 
dependent offspring and related descendants. This simply reflects 



the intensity of resource transfers embedded in human life history 
adaptations [41-48], set in large social networks with self- 
organizing spatio-temporally and functionally hierarchical struc- 
tures [49-53] (Fig. 4). In using his model to interpolate the 
population trajectory between 1 and 1995 CE — yielding a visually 
satisfying fit to intermediate historical estimates and modern 
census data — Cohen implicitiy supports the expectation that m 
does not vary across network scales, at least above a certain 
threshold in which socio-politically and economically complex 
cultural systems are constituted by geographically widespread 
metapopulations, with total census sizes greater than roughly 10 6 
or 10 7 . 

In the OC and CC models alike, the Marx parameter measures 
the individual's cost of biological reproduction in society, relative 
to investment needed to raise aggregate carrying capacity 
marginally. Both models are specified according to plausible 
general premises about how population change can drive niche 
construction, recursively altering resource availability in a way that 
may cause an elastic, disproportionately positive or negative 
response. Carrying capacity elasticity, c(t), is simultaneously a 
coefficient of population change, Pt+ht — Pt, and recursively 
determined by P t \ consequently, carrying capacity changes non- 
linearly over time, and the long-term coupled dynamics of P t and 
K t are highly sensitive to initial conditions. This logically raises the 
possibility that — in the context of the historical, roughly scale- 
invariant Marx parameter definining r soc relative to r max — late 
prehistoric and early historic biocultural evolutionary changes 
influencing the potential for economies of scale, L, could have first 
delayed impacts on carrying capacity and population, only to give 
way to modern accelerating, super-exponential growth. 

Methods 

Multiple historical estimates and recent United Nations census 
data provide a widely accepted reconstruction of long-term global 
population change. These data and references to the original 
sources are openly available on the United States Census Bureau 
website [24] and are shown in Table 1 and on Figs. 1 and 2. This 
section describes a method of sensitivity analysis of the fit of OC 
and CC model trajectories to the independently estimated 
historical data and modern census values. All calculations and 
analyses based on the methods were carried out in Microsoft 
Excel. 

Coupled Difference Equations 

As Cohen points out [22,23], it is useful to model coupled 
population and carrying capacity change with difference equa- 
tions, in order to explore complex dynamics across annual or 
generational intervals. Cohen writes the discrete-step difference 
equation for logistic growth: 

P,+A,-P t = rP t (K t -P t )At (3) 

In constructing the discrete-step difference equation for 
carrying-capacity change, we substitute rP t (K t — P t ) for 
(P t+At -P,)/At, so that K l+ht -K t = c,rP t (K t -P t )At. Following 
the definition of c(t) in equation (1), we write: 

K t+At -K, = r(L-aP h t + v )(K t -P t )At (4) 
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Figure 3. Northern Hemisphere temperature variation, 700 AD - 1 900 AD, juxtaposed with historical population estimates. Historical 
climate proxy data are from the supplementary materials in ref. [35]. Historical population data as in Fig. 1. 
doi:1 0.1 371 /journal.pone.01 05291 .g003 



From the equivalence — ^ — = c , the population elasticity 

of carrying capacity — that is, the responsiveness of carrying 
capacity, K, + ^, — K t , to a given change in population, 
P, + &,—P, — is simply c t [see equation (1)]. 

Population trajectories following the CC model may be 
generated from initial conditions by coupling equations (3) and 
(4). The OC model may be generated from initial conditions by 



calculating values of carrying capacity, population, and change in 
those variables by setting parameter a = 0, and coupling equation 
(3) with the CC model's equation (4). 

Sensitivity Analysis 

The historical estimates for global population in the year 1 CE 
range from 170 to 300 million. The United Nations and the 
United States Census Bureau use country-specific census counts 
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Figure 4. Typical preindustrial flow of transfers in human groups among all life history stages, emphasizing the importance of 
giving and receiving for prime adults engaged in family formation, investment in fertility, and offspring and descendant care. 

Transferred resources include consumable calories, material capital, social capital, technological know-how and environmental knowledge, and 
cultural competence. The preponderance of exchange within life history stages occurs among married adults and involves a combination of material 
and social capital, including reputation. In turn, married adults account for the bulk of transfers to other life stages. Modified after ref. [37]. 
doi:1 0.1 371 /journal.pone.01 05291 .g004 



and data on mortality and fertility rates to calculate that world 

population reached 7 billion in 2012. Sensitivity analysis can 

evaluate how well the OC and CC models interpolate global 

population trajectories between 1 and 2012 GE, in terms of their 

fit to independent census measurements within the 1-2012 CE 

interval. It can also measure the fit of the projected population in 

2100 CE to independently derived United Nations projections, 

which use trends in national rates of mortality and fertility. In 

order to trace how the OC and CC models behave, the sensitivity 

analysis considers six cases defined by a low and high Pj CE value 

(225 x 10 6 and 375 x 10 6 ), respectively, each combined with three 

separate r sm values (1.0xl0~ 13 , 1.5 xlO" 13 , and 2.0 x lO" 13 ). 

For each trajectory, K t , P t , and c t were calculated for annual 

intervals — from initial conditions at 1 CE to 2500 CE — recursively 

using equations (3) and (4). General fit of the resulting modeled 

time series was measured as the average distance, d, from the 

historical global population estimates and census data for the years 

p 

1750 to 2012 CE. Here, d, = \\og^\. For the estimated global 

Pi 

population intervals from 1750-1950 CE, P t is the arithmetic 
average of all measurements for a given year. For the UN 
population values from 1955-2012, P t is simply the mid-year 
global census estimate. The year 1750 CE was chosen as the 



beginning of the measurement of fit between modeled P t values 
and independently estimated P t values, because national censuses 
began to be made in a widening number of countries around that 
time [2,23]. Thus, the period from 1750 to the present includes 
increasingly precise and accurate population measurements. 

Preliminary analysis determined that values of m«lxl0 12 
(yielding values of r Jor «l x 10~ 13 ) generate population trajecto- 
ries more closely concordant with the historical estimates and 
census data than the smaller value (m « 1 x 10 8 ) used by Cohen. In 
fact, the OC model trajectory shown in Fig. 2 yields a distance 
value rfi75o-20l2C£ = 0.126; in comparison, best-fit trajectories 
based on the higher Marx parameter reduce ^1750-2012C£ by a 
factor of roughly three or four. In addition, values of L — the limit 
to economies of scale — that are larger than Cohen's illustrative 
value of 3.7 billion also appear to offer better-fit trajectories. 
Holding m constant, the complete annually resolved OC 
trajectories were calculated for values of L from 5-200 billion. 
For each OC trajectory calculated in each "L-scenario", the d- 
value was recorded, along with P isno ce, P2012 ce, P2100 ce, ar >d 
P2500 ce- Fig. 5 demonstrates how the value of L that minimizes d 
is sensitively dependent the scenario's initial conditions in the OC 
model. The OC trajectory that minimized d for each case 
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provided the value of L used to evaluate the CC model's fit to the 
historical estimates and census data. Holding L and r soc constant, 
combinations of a and b values were used to generate new 
population time series. The parameter a was varied between 
1 x 10~ 3 and 1 x 10~ 8 by whole orders of magnitude. For each 
level of a, the parameter b was varied between 0.1 and 2.5 in 



increments of 0.1. Finally, the value of b that contributed to 
minimizing d in the coarse analysis was used as the center of a 
finer examination of best fit across values of a between 0.1 and 2.5. 
Again, for each CC model trajectory, the <i-value was recorded, 
along with Pjsoo ce, P2012 ce, P2100 ce, an d P2500 ce- 



Table 2. Sensitivity analysis results for the Original Cohen (OC) model. 





P, CE 


*soc 


OC Model (no consumption costs for growth in 


carrying capacity) 






L 


£/ alI? (1750-2012 CE) 1 


P 1800 CE 


^2012 CE 


^2100 CE 


2.25 x10 s 


1.00x10"" 


6.1677x10 10 


0.048 


6.1026x10 s 


8.1964x10 9 


7.8536x10'° 




1.50x10 -13 


3.9684x10'° 


0.044 


6.3341 x10 8 


7.7096x10 9 


5.8960x10'° 




2.00x10"" 


2.9045x10'° 


0.042 


6.5758 x10 s 


7.6053x10 9 


4.9666x10'° 


3.75 x10 s 


1.00x10"" 


6.0711 x 10'° 


0.030 


8.6160 x10 s 


6.9992 x10 s 


4.4469x10'° 




1.50x10"" 


3.9112x10'° 


0.033 


8.9298 x10 s 


6.7589 x10 s 


3.6172x10'° 




2.00x10"" 


2.8640x10'° 


0.035 


9.1781 x10 s 


6.5800 x10 s 


3.0971 x10'° 



'The value in this column is the average distance, d t — as defined in the text — between the OC model population size and the estimated or measured population size for 
year t from 1750 to 2012. Only the estimated or measured years, shown in Table 1, were included in the calculation of rfi75o-20i2C£. 
doi:1 0.1 371 /journal.pone.01 05291 .t002 
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Results 

Tables 2 and 3 summarize the variation in the OC and CC 
models' fit to historical and census data, respectively. In each of 
the six OC cases examined, there is a value of L that markedly 
minimizes deviations from the recent historical estimates and 
census data for the interval 1750-2012 CE. From this perspective, 
the OC model does indeed provide good fit with historical 
population data, closely tracking modern super-exponential 
growth over the past two centuries. However, these best-fit OC 
population trajectories involve limits to economies of scale, L, on 
the order of 10's of billions of people. Although the best-fit OC 
model successfully interpolates the population trajectory from the 
1 CE mean estimate (225 million people) to the 2012 CE global 
human numbers, it no longer predicts the slowdown in population 
growth that has been confirmed by the past two decades' census 
and demographic rate data (see Fig. 1). 

In fact, the historically best-fit trajectories of the OC model 
project global super-exponential growth to continue unabated 
through the 21 s century. For example, the OC model with 
parameters set as in Fig. 6 yields a population in 2100 CE at 59 
billion. The Malthusian steady state is reached by 2500 CE, at 286 
billion people. When carrying capacity, K t , can grow for so long 
without resource depletion or other systemic constraints, the 
velocity of population only slows down after two or three doubling 
periods past L. Only at this stage do population growth and 
carrying capacity growth converge asymptotically toward zero. 
Needless to say, 286 billion exceeds any of the historical and 
modern estimates of global human population saturation that 
Cohen reviewed [23]. 

In contrast, although the fit of the CC model to historical 
estimates and census data is generally similar to that of the OC 
I o model, it further achieves remarkable agreement with independent 

demographic projections of declining population growth through 
the 2 1 st Century. Fig. 7 illustrates how finer variation in parameter 
values and initial population and carrying capacity conditions may 
be tuned to generate closer agreement between the CC model 
■| ° trajectory and the 21 st Century UN demographic projections. 

.g 2 Discussion: The Costs of Consumption in Human 

Population Systems 

The sensitivity analysis of the original Cohen model and the 
consumption-cost extended version supports the hypothesis that 
the super-exponential population growth of the 19 th and 20 th 
centuries was only proximately caused by modernity's organiza- 
tional, ideological, and technological changes. Although only one 
^ £ u line of evidence, the sensitivity analysis of the OC and CC models 

3> 3 -j= suggests that recent population growth did not involve a simple 

£ .i?o and sudden, technology-, ideology-, or human capital-dependent 

- .2 transition from Malthusian constraints. Those constraints would 

J C ti 

& > have remained. However, during the period of rapid global 

J ~° m population growth, the marginal gains in carrying capacity, K t , 

^ 15 oi °. yielded by economies of scale continued to exceed the costs of 

^ ^ ^ Si adding more mouths to feed (see Figs. 6 & 7). 

™ 2 c 5 This result potentially illuminates why Malthus did not 

I >. | jB anticipate the modern human population explosion. Malthus 

S g 8 q- incorrectly assumed that food supply was mainly dependent on 

area under cultivation. 19th and 20th Century aggregate food- 
supply increase outstripped population growth, not only due to 
rising agricultural yields, but also due to organizational, legal, 
transportation, processing, storage, educational, and ideological 
innovations — all of which were network-scale dependent and 
mutually interdependent [54]. Recent theoretical treatments of 
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YearCE 



Figure 6. Consumption-cost (CC) model trajectory for the case involving i- JW = 2.0 x 10 13 ; « = 5 x 10 u ; A = 2.04; P\ C e = 2.25 x 10 s ; and 
K iC E = 2.26 x 10 s . See Table 3 for additional details. 
doi:1 0.1 371 /journal.pone.01 05291 .g006 



modern demographic dynamics emphasize how material and 
energy-intensive intergenerational resource-transfer strategies have 
simultaneously favored rapid fertility declines and rising consump- 
tion rates [10-12,55,56], but the OC and CC models alike suggest 
that modernity's most marked biocultural trade-off — reallocating 
resources from investment in fertility to per capita transfers — has 
not significantly increased global aggregate potential for achieving 
economies of scale. Rather, substantially earlier biocultural 
evolutionary developments — likely involving previously unrecog- 
nized systemic constraints on how family formation practices and 
intergenerational transfer norms structure and are structured by 
wider economic coordination and competition over extraction, 
production, and distribution — would have set limits to growth in 
the potential for achieving economies of scale, represented in the 
OC and CC models by L. 

It must be remembered that these models, as presented, involve 
explicit simplifying assumptions about both the potential for 
economies of scale, L, and the intrinsic social-network-mediated 
rate of increase, r soc (mainly shaped by the Marx parameter, m). 
These factors are hypothesized to be essentially constant features 
of historical, socio-politically complex, geographically widely 
interconnected cultural systems and the populations that constitute 
them. The general coupled difference equation approach to 
population and carrying capacity change may prove helpful for 
studying regional — rather than supra-continental or global — 
demographic dynamics. Here, though, more detailed parameter- 
ized models will likely yield greater insight into demographic 
history variation among such relatively local contexts [57-61]. In 
this study the OC and CC models are scientifically relevant 
because they evoke plausible supra-regional contexts in which the 



parameters L and r soc are roughly constant. Here, it becomes clear 
that the OC and CC models are scientifically useful only insofar as 
they help us think through and refine relevant hypotheses about 
very complex interactions among human populations, the social 
networks they constitute, and the environments with which they 
extract, transform, and exchange matter and energy. From this 
perspective, inspection of the results for the CC model suggests 
that we may be able to trace possible cultural evolutionary 
transitions in levels of L, r soc , and K(t). For example, comparison 
of the six cases analyzed in Tables 2 and 3 shows that — for the 
OC and CC models alike — the higher value of Pj ce (375 million) 
yields the better fit to historical population estimates. This initial 
population value is also consistent with a more complex 
demographic history, in which population growth from the first 
millennium CE through the Medieval Warm Period (see Fig. 3) 
may have involved a substantial increase in the ideological, 
institutional, and technological foundations for potential econo- 
mies of scale, L. Although the Marx parameter, in — that is, the 
average cost of forming a family, investing in fertility, and 
obtaining and holding transfers to offspring — would have 
remained very high compared to r nulx , it may also have declined 
significandy during the apparent Medieval Warm Period era of 
demographic growth. 

In all instances, key political-economic changes — involving 
organizational, ideological, and technological innovations — would 
have shaped the potential for economies of scale, supporting very 
high population elasticity of carrying capacity around ca. 1500 
CE. Indeed, in the closely fit CC model trajectories shown in 
Figs. 6 & 7, carrying capacity, K t , begins increasing faster than 
population, P t , around 1500 CE. Yet, the extraordinary fit of the 
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Year CE 

Figure 7. Consumption-cost (CC) model trajectory tuned by inspection to fit P2012 c£= 6 - 9 " billion. In this case, r„=2.0x 10 13 ; 
L = 6.0779x 10 10 ; a=lxl0-> 4 ; b = 2.54; P XCE = 3.75000 x 10 s ; and K XC e = 3.75966 x 10 s . This case results in 5i7 5 o_2oi2 = 0.031 and 
Pi\mce = 1.584 x 10 10 , demonstrating that — according to the CC model — initial conditions established between 1 and 1500 CE can largely account 
for the historical pattern of modern super-exponential growth and the projected trend of population deceleration through 2100 CE. 
doi:1 0.1 371 /journal.pone.01 05291 .g007 



CC model trajectory to historical estimates, census data, and 21 st 
Century demographic projections suggests that critical political- 
economic developments were in place much earlier. Whatever the 
exact nature of these prior changes in organization, ideology, and 
technology, they were initially adopted in biocultural evolutionary 
environments in which population growth was limited over large 
continental or supra-continental land areas (that is, >ca. 10 1 - 
10 8 km 2 ), more often than not subject to growth-limiting 
Malthusian constraints over generational to centennial timescales. 
The CC model is highly sensitive to initial conditions, and the 
best-fit trajectories all involved initial carrying-capacity levels (K 0 ) 
only slightly greater than actual initial population levels {Pq). In 
other words, it appears that even with beginning conditions 
incorporating very high limits on economies of scale, L, those same 
initial conditions would have also had the global human 
population just barely under the prevailing carrying capacity — 
that is, very close to immediate Malthusian limits on growth, 
Po»Kq. This would establish long-term suppression of population 
rise, while contributing to a nonlinear, positively elastic carrying- 
capacity response. Significant cultural structures determining 
population-dependent limits on achieving economies of scale 
would have emerged when continental-scale populations were 
largely stationary over a period of many generations or centuries. 
The historical estimates suggest that global population growth was 
stagnant during the entire first millennium CE, although 
archaeological and historical evidence documents comparably 
dynamic variation in migration and settlement patterns, on the 
one hand, and political economy, on the other [38,62]. Here, the 



OC and CC models help to direct our attention toward a new 
hypothesis. Developments in cultural system complexity — emerg- 
ing in different parts of the world in the centuries and millennia 
prior to 1500 CE — were more systemically important in setting the 
stage for the recent population explosion than was any specific 
modern technology or ideology introduced during or after the 18 th 
Century. 

How then might unambiguously pre-industrial and pre-modern 
biocultural evolutionary processes have generated an enormous 
potential systemic capacity for achieving what turned out to be 
19 th and 20 th Century industrial and information economies of 
scale? The answer arguably lies in dynamics of competition, 
power, and inequality [10,38,63,64]. Let us assume a prehistoric 
set of initial conditions — emerging during the Holocene, between 
the origins of agriculture and the establishment of states and urban 
settlement systems — in which substantial within-population het- 
erogeneity in political power and biological well-being have 
become institutionalized, as continental-scale meta-populations 
approached a stationary steady-state. Under such conditions, 
Malthusian regulation would have largely prevailed. However, 
within populations an elite segment would have been able to 
exploit differential access to material resources, information, and 
media of communication, in order to mobilize labor on risky or 
expensive organizational and technological innovations. Reflecting 
the extant variation in access to information, ideological interests, 
and power, a large number of such innovations — although they 
would have had organizational and technological forms — would 
have actually had ideological, political, and military functions. 
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Figure 8. Population elasticity of carrying capacity c t (see equation 1) and % annual population change over time, for the case 

r soc = 1.0 x 10 °; £ = 6.0711 x 10 10 ; « = 2x 10 14 ; 6=2.20; P 1C e = 3.75 x 10 8 ; and K 1C B = 3.76 x 10 s . As population change reaches an equilibrium 
steady state, population elasticity of carrying capacity is substantially negative, revealing the equilibrium to be weak. 
doi:1 0.1 371 /journal.pone.01 05291 .g008 



In this setting, competition for power and material resources 
could lead to an increase in the capacity for achieving economies of 
scale, L, even as carrying capacity, K t , persisted at or close to 
prevailing population levels. This would occur as competition, 
exploitation, domination and resistance processes favored increas- 
ing political complexity and multiple shifting paths of economic 
connectedness, shaping the emergence of heterarchically and 
hierarchically related institutions [65,66]. These socially structured 
and structuring institutions and activities would effectively 
consume biological well-being — that is, fitness in the context of 
transfer-intense, extended life history strategies (cf. ref. [20]). This 
would initially limit population growth. Here, political and 
economic factors would alter the scale and distribution of variation 
in biological well-being, while depressing carrying capacity, K t . 
Intense competition among institutions, actors, and shifting 
alliances would have simultaneously depended on and limited 
growth in material and labor resources. Cultural evolution may 
have increased the physical inputs, flow rates, and outputs (in 
terms of economic production and fertility), but the net 
demographic effect would have been near zero. 

The organizational and ideological systems emerging from such 
cultural selection would have, in turn, increased the logistical, 
ideological, and technological limits on achieving economies of 
scale, L, without significantly raising population carrying capacity, 
K t . Only later, these cultural evolutionary processes — including 
diversification and competition among ideologies; economic 



system complexity and resilience; and proliferation of political, 
religious, military, and economic institutions — would have sup- 
ported the development of positive feedback between population 
and carrying capacity. 

Conclusion 

In comparative and evolutionary perspective, recent human 
super-exponential growth seems unlikely — or at least unfamiliar 
under standard theoretical models. Darwin's early insights about 
natural selection were crucially influenced by his reading of 
Malthus. In his Notebooks on the Transmutation of Species [67,68], 
Darwin wrote (sic): 

Even a few years plenty, makes population in Men increase & 
an ordinary crop causes a dearth, take Europe on an average 
every species must have same number killed year with year by 
hawks, by cold &c. — even one species of hawk decreasing in 
number must affect instantaneously all the rest. — The final 
cause of all this wedging, must be to sort out proper structure, 
& adapt it to changes. — to do that for form, which Malthus 
shows is the final effect (by means however of volition) of this 
populousness on the energy of man. One may say there is a 
force like a hundred thousand wedges trying force every kind 
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of adapted structure into the gaps in the oeconomy of nature, 
or rather forming gaps by thrusting out weaker ones. 

Here, Darwin took clear note of Malthus's insight that 
"geometric growth" yields startlingly short, decadal-scale human 
population doubling times. This underpinned Darwin's elegant 
argument that normal intra-decadal-scale environmental fluctua- 
tions should drive recurrent episodes of density-dependent 
competition within animal populations. This should limit long- 
term population increase but maintain conditions for what he 
came to call natural selection [69] . In this light, the rapid human 
demographic growth of the past 200 years is quite simply a 
remarkable population biological phenomenon. The consump- 
tion-cost model presented in this article formalizes plausible 
conditions in which an apparent Malthusian trap — where steady- 
state population regulation is theoretically expected to inhibit 
political-economic risk-taking, thereby limiting the adoption of 
technological or organizational innovations [60,70,71] — actually 
belied a very different biocultural evolutionary situation. It was not 
a fixed environmental carrying capacity that held historical pre- 
industrial populations in check. Carrying capacity and population 
growth alike were instead temporarily limited by intense political 
competition, economic dynamism, and change in the hierarchical 
scale and heterarchical institutional diversity of the overall cultural 
system. Yet, the evolution of joint political and economic 
institutional complexity increased the potential for social-net- 
work-dependent carrying capacity growth. This sets up conditions 
for long-term population change involving a delayed, yet 
extremely rapid acceleration in population increase. 

The consumption-cost model's success as an interpolation 
function, for the period 1-2012 GE, thus builds on earlier insights 
from the original Cohen model [22,23]. The best-fit model 
population trajectories (see Table 2; Figs. 6 & 7) support the 
hypothesis that both modern super-exponential population 
increase and the incipient, likely rapid population deceleration — 
which we are currendy beginning to face — are part of a long-term, 
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